Quantum Statistical Mechanics for Nonextensive Systems II 



E. K. Lenzi^, R. S. Mendes^ and A. K. Rajagopal^ 
^Centra Brasileiro de Pesquisas Fisicas, R. Dr. Xavier Sigaud 150, 
22290-180 Rio de Janeiro-RJ, Brazil 
Departamento de Fisica, Universidade Estadual de Maringd, 
Av. Colombo 5790, 87020-900 Maringd-PR, Brazil 
^Naval Research Laboratory, Washington DC 20375-5320 
(February 1, 2008) 

Abstract 



In this paper, the Green function theory of quantum many-particle sys- 
tems recently presented is reworked within the framework of nonextensive 
statistical mechanics with a new normalized g-expectation values. This re- 
formulation introduces a renormalization of temperature of the earlier theory 
and a self-consistency condition. The linear response theory is also presented, 
along with its two-particle Green function version. Finally, a Boltzmann 

transport-like equation is also developed here. 
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I. INTRODUCTION 



In two recent papers [|T],0 a Green function theory of nonextensive many-particle systems 
was developed. This work is here reformulated in terms of a new version ||^ of the Tsallis 
ensemble involving normalized g-expectation values in the generalized context. This new 
version obtained in this way has the following desirable properties which the early versions 
lacked invariance of the ensemble with respect to uniform translation of the energy 

spectrum, the expectation value of a c-number is the same as that c-number and finally, 
the preservation of the formal additive structure of an operator associated with two subsys- 
tems (e.g., energy). All these features entail the appearance of renormalized temperature 
parameter replacing that in the previous formalism along with a self-consistent condition. 
The purpose this paper is to revisit the general Green function theory presented in |l|J^ in 
terms of the new framework. In Sec. II the one-particle Green function theory for both the 
Bose and Fermi systems along with explicit results for free-particle systems is given. These 
have implications to experiments as described before In Sec. Ill a generalization of 

the linear response theory ^ in the new framework is given and it is re-expressed in some 
special cases in terms of the two-particle Green-function. This is related to scattering cross 
section as before In Sec. IV a Boltzmann transport-like equation is, for the first time, 

derived, for arbitray q. The final Sec. V contains conclusions of this work. 



II. THE GREEN FUNCTION THEORY 



The Tsallis entropy Sq = {1 — Trp'^)/(g — 1), Trp = 1, p the system density ma- 
trix is maximized subject to normalized g-expectations values {H)q=Uq = TiHp'^/Tip^ , 
{N) g=Ng=TTN p'^ /Tt p'^ . A detailed discussion of this in relation to the earlier formulations 
may be found in [Q. To distinguish from our earlier work, a tilda over the corresponding 
expressions in the foregoing sequel is used. The new density matrix determined in this way 
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is then 



P = expg (-X(/?,/i)j /Zq , 



Zq = Tr 



(1) 



where 



1/(1-^/) 



exp^ {-X0,fi)) = (l - (1 - 



(2) 



with P = P/cq, and Cg = Trp"^. Here /3 is the normahzed temperature. /3 and are the 
Lagrange multiphers as before associated with the normahzed g-expectation values given 
above. The constant Cg is a function of all these quantities. There are two self-consistency 
conditions that follow from Eq.(P, because 



from which it follows that 



q)P((H-Ug)-fiiN-Ng)) 



c = Z^~'^ 



(3) 



(4) 



By using the expression for p in Eq.(|Ij) in the definition of Cg, another relation is derived: 



Cg = Tr 



expg 

Using Eq.(|D in Eq.(^ an equivalent expression is obtained for Zg-. 

Zg = Tr [expg /i)^^ 

Noting that exp^ (^—X{(3, jj,)^ = 
Tr' 

ogy. The expression for a general normalized g-expectation value is then 



(5) 



(6) 



exp^ X(/3, /i)^ |l — (1 — and since 

exp [— X(/3,/i)) '^X = by definition, these two expression for Zg become a tautol- 



{A)q = Ag 



TiA 



Tip 



exp^ 



Tr{i 


expg 
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Tr " 


expq ( 


-X0,p))' 
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(7) 



Here Eqs.(^^) were used. From Eq.(0), Uq, Nq are deduced which along with Eqs.([|,|,mD 
form the required self-consistent relations. In this form for the expectation value, only the 
"connected" diagrams appear and thus explain the results found 0,^] contract in to that 
obtained in |^ which did not used normalized g-expectation values. The one-particle Green 
function is then redefined as in with the new normalized g-average and the density matrix 
given above: 



G(^)(l,l';A/i) = -(T(vI/(l)M/t(l'))) 

iTr [3'T(^(l)^t(i')) 



9 ' 

(8) 



Using the contour integral representation employed in 

fei-f / du exp(-« h){-u)-^ = -1- (9) 
2n Jc I [z) 

with b > and Re 2; > 0, where the contour C starts from +00 on the real axis, encircles the 
origin once counterclockwise and returns to +00, Eq.® is expressed in terms of the Green 
function with q = 1: 

G('')(l,l';/3,/x)= / du^'\u)Z,{-u{l-q)P,^i)G^'\l,l';-u{l-q)P,^i) , (10) 
where 

kf\u) = .^Mi^exp(-.)(-.)Va-.) 

/TT Ziq 

1= / dwirf . (11) 

Note that all these definitions are consistent with Eq.(0) defined above. Note also that 
these expressions have similar appearance as those found in When q = 1, the constant 
terms involving, Ug, Ng factor out so that the one-particle Green function in integrand of 
Eq.(p!oD is the conventional extensive one: Zi appearing there however is the corresponding 
conventional partition function multiplied by e~'^^^^~'^''^^''~^^''\ From Eq.(|I]) an equivalent 
expression is found for Zg\ 



Z,= [ dnir«HZi(-n(l-g)/3,/i) 



ic 

i 

~27r 
Z, 



iV f^") 

^(1) ^ _kiziZexp(-n)(-M)-(2-9)/(i-a) 
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m(1 - q) 



kf\u) . (12) 



As in 0, correlation functions are introduced. Following KB ||T0|, we introduce correlation 
functions 

Gi«)(ll';A/i) = -(vi/(l)vl/t(i'))^, 

G^^\n'-(3,f,) = j{¥{l')^{l)),. (13) 

The notation > and < is intended to exhibit the feature that G^'^\l, 1'; /5, jj) = Gy\l, V; (3, jj) 
for ti>tv and G'('?)(l, 1'; /5, /i) = G^l\l,l'] (3, ^i) for h < ty. Using (D, we may similarly 
express G>'' and in terms of the corresponding grand canonical correlation functions. 
The spectral weight function in frequency space by taking the Fourier transform with re- 
spect to time differences, yl(ri, ri/; cj), introduced in KB reflects only the properties of the 
Hamiltonian H . The average occupation number in the grand canonical ensemble of a mode 
with energy oj, f{uj,(3) = (exp(/5(ti; — fi)) =F 1) ""^j takes account of the basic permutation 
symmetry of the system. We can thus express G'>'' and in terms of these in the following 
way: 

iG'>^(fi,ri/;cj;/3,/i) = j^(\ukf\u)Zi{-u{l - q)j3,n) iG>\fi,fr;uj; -u{l - q)j3,n) 
= I dwKf (w)(l±/(cu,-w(l-g)/3,/i))A(fi,n,;cu)Zi(-w(l-g)/3,/i) , (14) 



iG^l\f,, fy- u- P,fi)= duKj;'\u)Z^{-u{l - q)(3, fi) ^G^^\r^, r,,; -u{l - g)/3, /i) 

J c 

= ±f duKjj'\u)f{u;,~u{l-q)P,fi)A{n,r,,-,u;)Z^{-u{l-q)P,fi) . (15) 
J c 

Thus the spectral weight function is found to be 

A{ri,fi,]uj) = i (Gj^(fi,fi.;a;;/3,/i) - (^^^^(fi, fi/; /3, /x)) 

= / duK!;'\u)A{n,nr,uj)Z,{-u{l-q)^,fi) . (16) 



From this we deduce an important sum rule 



°o duj 



(17) 



This is just an expression of the equal time CCR of the particle fields. The spectral function 
weight is a property of the given systems not dependent on the ensemble and the results 
have similar appearance as in 0] with the modifications noted here. Using the Fourier 
representations of the step functions involved in the time-ordered Green function in Eq.(§), 
an important result, generalizing the result found in M is deduced: 



^ c 



oo 2tt 



oo 



XI —A{fi,fir,uj'){P 

-oo ZTT Vu; — W 



- 1 —i7T6{uJ — Uj') 



tanh ( - - lA 



;i8) 



Note that in the view of Eq.(|T6|), only the delta function term in this equation depends 
explicity of q as in Ref . [§] . So far we have discussed the one-particle properties. The above 
development is similarly extended to generalize the many-particle g-Green functions. Using 
the same notations as in KB, we have in general. 



G(^)(12...n,l'2^..n';/3,/i) = 4w(l)vI/(2)...vl/H^l/t(l0^1/t(20...v&n^'))), 
dukf\u)Z^{-u{l - q)(3,fi) G«(12...n, l'2'...n'; -u{l - q)P,fi) . 



(19) 

For a many-particle system with a Hamiltonian containing one-body potential Vi(ri) and in- 
stantaneous two-body interaction potential V2{ri,r2), which is symmetric under interchange 
of 1 and 2: 

V^^(r,t) ■ V^(r,t) 



H 



df 



2m 



+ / dr \/i(r)*^(f,t)^(f,t)) 



+ J drdf ' ^■f(f,t)^'^(r",t)\/2(^,r ')*(r",t)*(f,t) . 
We can express Ug in terms of the Green function following KB 



(20) 



U, = ^ Idf 



' d__d_ 
dt'df 



m 



Vi{r)~Vi{r') 



G'^£{f,t;f',t'] 



f'=rt'=t 



(21) 



For the free particle case Vi(r) = as in Eqs.(p3D, (]2^) we have 



y Jc " ^ ' J-oo2ti J {2ti)d \ 2 I (e-«(i-9)/3(c^-M) ^ 1) ^ ^ ' 



For a uniform system, we can take Fourier transforns with respect to fi —ry in Eq. (p!5D and 
express the one-particle momentum distribution function (iV(p ))q in terms of the spectral 
weight function of the N-particle system. 



Similarly the one-particle frequency distribution function {N{u!))q is given by 

f d^p 



(23) 



N{uj)g = ±V I dukf\u) 



(24) 



(^g-«{l-g)/3(a;-f.) ^ ;L) ^ {27t)^' 

Here V is the volume of the D dimensional space in which the particles reside. The chemical 
potential is determined by the expression for the g-mean value of the total number operator 
N , 



± / duirf)( 



u 



^ duj r d^p Zi{-u{l- q)(3,n) 



A{p- uj) 



(25) 



V Jc ^ '' 'J-oo27tJ (27r)^ (e-"(i-9)/5{'^-^) ^ 1)' 
Equation (p3D will now be used to derive the one-particle g-momentum distribution func- 
tion for a free fermi gas, where A{p, uj) = 2116 {u — e^), being the one-particle energy. The 
known partition function for the free fermi gas in the Gibbs theory need in this expression 
is approximated in a simple form 



^ \ 3/2 00 ^f3fj.l 



2np) g 

and hence also a further approximation resembling the fugacity expansion of the classical 
Maxwell gas. 



V 



m \ 
2^ 



3/2 



(26) 



00 -| 

^i(/5,/^) = E7r 

1=0 



V 



m 



27rp 



3/2 



(27) 



and l/(e^ + 1) = Z^jLil^l) e for |e < 1, we have finally the result. 



V 



00 



1=0 



V 



m 



271 f3, 



3/2 



j= 



31 



{1 + (1 - g)/3 {Ug -ixNg + ixl-J {ep - /i))} 



g/(l-g)+3//2 



[28) 



An approximate expression for the g-mean momentum distribution function in the first 
version of the Tsalhs formulation [Q] (where ordinary g-mean values were used) has been 
used in the literature |TT|. We can deduce such an expression as in |TT[ from the above by 



making a further drastic approximation by taking / = 0,Zq = 1, dropping Uq, Nq and (3 = (3 
to be in conformity with the first version of the Tsallis formalism [Q, we obtain the result, 



^^f;(-l)^+i{l-(l-g)/3J(e,-;.)r/(^-) . (29) 
^ j=i 

A futher approximation corresponding to the high temperature scheme where /3 (e^ — /i) << 
1, we have 

(30) 



In order to assess the accuracy of the approximation in Eq. ([30| ) in relation to the more exact 
result in Eg. (pQ]) , we compare the respective second order terms of the series, {qJ{qJ — 



[1 — q))/2) [13 (ep — /i))^ and (g(2g — l)/2) (J/5 (ep — /i))^, and we find the error involved is 



of order (1 — g). This result was essentially obtained in |[TT]] via a different approach to this 



problem. Similar calculations for the Bose and the Maxwell gases may be made on the same 
lines as above yielding corresponding expressions for the respective momentum distribution 
functions. 



III. LINEAR RESPONSE FUNCTIONS 



The linear response is now calculated, following [Q, using this new prescription: 



(A5(t)), = Tr P{t)B 



Tr 



P{H;qJ)B 



(31) 



where P{t) = ffl{t) /Tiffl{t). With this new definitions, all averages are calculated with 
proper unit normalization. The time- dependent P-operator in Eq.(^ID obeys the usual 
equation of motion with time- dependent Hamiltonian H — AX{t) 
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i^^P(t) = [H,P{t)] - [A,P{t)] X{t) 



(32) 



with the initial condition P{t = — oo) = P{H; q,P) given by Eq.(|l]). Taking the trace over 
both sides, we observe that TrP(t) = 1 for all times. The solution to this linear order in 
X(t) is the found the standard procedures: 

1 /•* 



P{t)^P{H;q,(3)- 



ih 



dte-*(*-*')^/" 



A,P{H;q,(3) 



^i{t-t')H/n 



(33) 



Thus 



(AS), 



d0M-t')X{t') 



ih 



Tr 



A,P{H;q,f3) 



B{t) 



ih 



-Tr 



A,B{t) P{H-q,P) 



(34) 



It appears therefore all the results and conclusions of [^] hold with the new definition of 
mean values with the renormalization of the temperature and the self-consistency condition 
in Eq.(|,|). 

We now turn our attention to rewriting the dynamic response and the scattering cross 
section in the g-formalism in terms of the integrals over the usual ones as was done above. 
Reformulating the result obtained in Ref. 0, the dynamic linear response of a quantity B 
to an external probe that generates A in the g-formalism is 



Xif]i(^,/3,/i) = hm /" dt 



"^^*-^*-ftl(t,/5,/x) 



(35) 



This in terms of the integral representation, is: 



xa(o;,/3,/i)= / dn^f («)Zi(-n(l-g)/3,/.)xa(^,-«(l-g)/3,/^) 
J c 



(36) 



where Xba{^^ —u{1 — q)f3,fi) is the usual Kubo dynamical response function evaluated now 
at a temperature —u{l — q)(3. Following Ref. in the context above and rewriting the 
g-averages of the anticommutator and commutator expressions, we have: 
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P{H, N; q, (3, /x) [A{Q)B{t) + B{t)A{Q) 



(37) 



$gl(t,/5)=lim/dt'e-^*'Tr 



P(i7,iV;g,/5,/i)[A(0),5(t) 

<j C 



(38) 



The fluctuation-dissipation theorem due to Kubo ||T^ for the extensive case, (g = 1) is 



(1) 



with 



7-. / \ ^ 1 / 

Efi{(^) = -coth I — 



(39) 



(40) 



Here we obtain 



X coth (-u{l - q)p) $gi(^, -u{l - q)P,fi)) 



(41) 



We now relate the scattering function defined for example in, Lovesey , in the g-formalism 



as 



dtexp(-zcjt)(i"f(o)i(t))W 



(42) 



2tT j-oo 

where A is the operator which affects the change in the states of the system in a scattering 
process. Here the superscript (c) denotes canonical ensemble instead of the grand canonical 
ensemble used earlier. This is equivalent formally to setting /i = in the earlier development. 
Then, using our transformation, we express this scattering function in terms of the usual 
q = 1 scattering function 



, io, —u{l — q)f3) . 



c 



(43) 



From 1^, by taking B = A\ we have that the imaginary part of the g-susceptibility, 
X^t^(t^,/5) can be expressed in terms of the q = 1 scattering function 
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X (l-exp(-M(l-g)/3)5W(A?,cu,-M(l-g)/3) . 



X 



(44) 



We have thus expressed the g-scattering function as well as the imaginary part of the asso- 
ciated g-susceptibility in terms of the parametric integrals over a kernel multiplied by the 
usual scattering function which now depends on this parameter as displayed above. It is 
possible to obtain other properties as in Ref. to this new version of the Tsallis statistics. 



There exists a class of disturbances which are not conveniently described by the usual 
equilibrium Green's function |]10[. For example, the disturbance produced by the externally 
applied force field, -F(r, t) = VU{f, t). And many interesting physical phenomena appear as 
response of systems to external disturbances of this kind, for example in an ordinary gas, a 
slowly varying U{f, t) produces sound waves. Thus, we expect that this feature will appear 
in the generalized g-ensemble theory as in g = 1. It is interesting therefore to extend the 
usual Boltzmann equation to the generalized case to understand this in the nonextensive 
context. In our formalism, this disturbance (force) may be represented by. 



where n{r,t) = '^^{r,t)'^u{1",t)■ Let us to obtain the colisionless Boltzmann equation from 
the Green functions given above using the Hartree approximation in the same spirit as KB. 
By using the Heisenberg notation we have that the g-expectation average is given by 



IV. BOLTZMANN TRANSPORT-LIKE EQUATION 




(45) 



Ei pUhto\Xu{R, t)\i,to) 



(46) 



where 



Xu{R,t) = V{t)-'X{R,t)V{t) 



(47) 



with 
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V{t) = T |exp 
From the Eq.(|l^) we have that, 



-t J dt' I d-'r'n{r',t')U{r',t') 



(48) 



1 



(49) 



G'(«)(12,l'2';f/) 



(T(*f;(l)^'^(2)vl/t^(2)vl/t^(l))) 



(50) 



In terms of these Green's functions we may describe the response of a system, initially in 
thermodynamic equilibrium, to an applied disturbance U{f,t). In the same way as above, 
the average density at point R, T is 



{n{R,T))u,,= {^lj{R,T)^u{'R,T))g 
= ±iGj^(R,T, R,T; U) 



(51) 



and the current at the same point is 

■ V - V r 



{J{K,T))u,, 



2mi 



±iG'^^\R, T, K',T;U) 



(52) 



R=R' 



where the conservation laws for the number of particles, the energy and the momentum are 
preserved here as well as in the usual case, and we can use them for derivation of the sound 
propagation. Here we define /^(p, R, T) (with r = ri — fii, R = l/2{ri + ry) t = ti — ty 
and T = l/2(ti + ti/)) as. 



/,(p,R,T) = J ^Gj)(p,c.,R,T; U) 



(53) 



generalizing the original definition proposed by Wigner. As in the usual case /g(p, R, T) 
leads to the generalized g-particle density 



(2vr) 

and the generalized g-particle current 



(54) 
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{AR,T))u,, = J^s^fM^^T) . (55) 

This identification above of tlie distribution function /^(p, R, T) will enable us to see the re- 
lationship between to Green's functions, transport equations, and the colisionless Boltzmann 
equation. Thus, to obtain the colisionless Boltzmann equation we use 

' 9 V? ^ 



± [ [ dt2dr2vin - r2)5{h - t2)G(^)(12, 1'2+; U) 



(56) 



and 



G('?)(12,l'2';f/) = G('')(l,l';f/)G'('')(2,2';t/) 



(57) 



After some simplification, we obtain 



^g^^+^-u,„il)\&'Hl^'■,u) = s^l-y) 



where 



(58) 



t/e//(R, T) = f/(R, T)±tJ dR'v{R - R')g1?^(R', T; R', T) . (59) 
By taking the difference of two equations in the variables 1 and 1', we find 

Considering now tii = tf = T and expressing Eq. (|60D in terms of, r = ri — ri/, R = 
1/2 (ri + ri') and using Eq.(p3|) we have that 



l7e//(R+^,T)-f/e//(R-^,T) 



^ / g ^ Vr ■ Vr 1 

I dT im i 



X 



(61) 



Now let us suppose that Ueff(R., T) varies slowly in R, then we obtain after some calculation 
the coUisionless Bolztmann equation 
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d p 



Vr- VRf/e//(R,T)- Vp 



/g(p,R,T) = (62) 



where 

U,ff{K,T) = U{R,T) + JdK'v{R-R') J i!^/^(p', R', T) . (63) 



V. SUMMARY AND CONCLUSIONS 

In this paper we have developed in detail the Green function theory for nonextensive 
systems based on the g-ensemble of Tsallis. By means of a contour representation, Eq.(P), we 
have made this theory resemble the usual one for extensive systems given by L. P. Kadanoff 
and G. Baym ||10[, for example, even though in actual practice, the results are very different, 
as exemplified by the representative results given in Sec. II for a variety of situations. Before 
this development, thermodynamic quantities for model systems were computed in the Tsallis 
ensemble as for example in Ref . P,p!^ . With the present work, we believe that the theory of 



many particle systems for the Tsallis ensemble has been considerably extended and placed 
on par with the conventional theory based on the Gibbsian ensemble, in that we have been 
able to compute response functions besides the thermodynamic quantities. The case of q 
different from unity is expected to apply for long-range interacting Hamiltonian systems 
0], among others. From a formal point of view, noninteracting and short-range interacting 
systems are mathematically well posed problems only for q < 1. In conclusion, we have 
here developed the formalism associated with Tsallis statistics for describing nonextensive 
many-particle systems by a suitable generalization of the corresponding Green function 
techniques, so commonly employed in such studies for extensive systems. As with the usual 
Green function theory, which has been traditionally successful in explaining experimental 
observations, the present work enables us to analyze future possible experimental work on 
nonextensive systems. We may add that other forms of the nonextensive entropy such 
as S^""^ = Sg/cg, deduced from considerations of form invariance of the statement of 
maximum entropy principle and the metric structure in quantum density matrix theory, 
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many also be expressed in contour integral form. The work presented here thus admits of 
applications elsewhere. 
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